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ABSTRACT 

Wide separation leasing statistics offer information about the density profile and abun¬ 
dance of dark halos. Recently a possible discovery of six quasar pairs, which may be 
leased multiple images, was reported by Miller et al. (2003). These pairs are selected 
from a catalog of the 2dF quasar (QSO) survey comprising 22163 quasars. We calcu¬ 
late expected numbers of leased quasars taking account of the redshift and magnitude 
distributions of the quasar catalog. Given some of the six pairs are genuine leased 
systems, we put interesting constraints on the inner slope of dark halos, Oq, and 
We show that the detection of even one lens with separation > 30" is marginally con¬ 
sistent with models that have cuspy inner density profile, p oc and very large 

csj CTg > 1.2 for Oq = 0-3. To reconcile with constraints from X-ray clusters or cosmic 
shear, much lower and much higher erg are needed, although such high erg seems 
too extreme. The shallower inner density profile p (x r~^ is hardly acceptable. In par¬ 
ticular, the expected number of lenses with separation > 200" is too small to explain 
the discovery of such anomalously wide separation lens systems. These results imply 
that we miss some important systematic effects, there is a problem in the cold dark 
matter scenario, or none of these six quasar pairs is likely to be lensed images. 

Key words: cosmology: theory — dark matter — galaxies: clusters: general — 
gravitational lensing 


1 INTRODUCTION 


The Cold Dark Matter (CDM) scenario predicts the ex¬ 
istence of cuspy dark halos, and thus is expected to pro¬ 
duce significant numbers of wide separation lenses {6 > 6"). 
Statistics of such wide separation lenses are known to be a 


dark halos llMaoz et all ll 997: 

Wvithe. Turner. & Soereell 


Sarbu. Rusim fc Mai koo j 

iTakahashi & Chiball200ll: iLi & Ostrikedl2002HOffurill2002J). 


Although a number of radio surveys has tried, they 
could not find wide separation lensed quasars (e.g., 
Phillips, B rowne , fc Wilki nsorJl2001al : IPhillios et ahlEoPlbl: 
Ofek et alJl200ll.l2nn2ll .' The lack of wide separation lenses, 
however, does not conflict with the CDM scenario because 
the expected lensing rate is significantly smaller than that 
of small separ ation lensing (0 ^ 1 "). 

Recently. iMiller et alJ ll2003^ reported that they found 
six quasar pairs which may be lensed multiple images in the 
2dF quasar (QSO) catalog comprising 22163 quasars, al¬ 
though the number of lensed quasars may be even larger 


than this because they have done follow-up for only 11 
quasar pairs among 38 quasars selected as possible lens can¬ 
didates. These systems were identified to be lensed images 
from the detailed comparison of quasar spectra. The sep¬ 
arations of all these quasar pairs are larger than 30", and 
the separations of some pairs reach even ~ 200". If this sur¬ 
prising result is true, it offers a lot of information about 
dark halos. In this Letter, we calculate the expected num¬ 
ber of arcminute-separation lenses {9 > 30") in the 2dF 
QSO survey. We use the realistic density profile predicted 
in the CDM scenario, and take account of the redshift and 
magnitude distribution of the quasar catalog. Therefore our 
results can be directly compared with the observation. We 
show that the existence of such wide separation lenses in the 
2dF QSO catalog is marginally consistent with the “con¬ 
cordance” cosmology if the value of only erg is very large, 
erg > 1.2 for Do = 0.3. This means that we can strongly con- 
str ain dark halo prope rties if some of quasar pairs reported 
bv iMiller et alJ (l20n3h are truly gravitational lens systems. 
Throughout the paper, we assume a fiat universe Dq-I-Ao = 1 
and h = Ro/(100kms“^Mpc“^) = 0.7. 
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2 CALCULATION OF PROBABILITY 
DISTRIBUTIONS 


The lensing probability distribution at wide separation re- 
flects the properties of dark halos, rather than ga laxies 
iNakamura fc SiitolIl997l : iKeetorJIl998l : lOguril 1200A . The 
halo density profiles predicted by recent N-body simulations 
may be parameterize d as a one-parameter family, t he gen¬ 
eralized NFW profile lIZhaolllQ^ l.ling fc Sutoll200(ll : 


p{r) 


pcTit(^z')Sc{z') 

(r/r.r 


( 1 ) 


where pcrit(2) is a critical density at 2:. While the correct 
value of a is still unclear, the existence of a cusp with 
1 < a < 1.5 has been established in recent N-body sim¬ 
ulations IINayar ro.Frenk^ fc Whit^|l99d [Moore et aLlll999l: 
l.ling fc SutoIbnOlirlDiiaiSiigrir^^chiolDoi^ . For definite¬ 
ness, in this paper we consider two cases, q = 1 and 
a = 1.5, which cover the range of the CDM predictions. 
The scale radius Ts is related to the concentration parame¬ 
ter Cvii{M, z) = rvii{M, z)/rs{M, z). Then the characteristic 
density Sc{z) is given in terms of the conce ntration param¬ 
eter (see e.g.. lOguri. Taruva. fc Sutoll200lll . The lens equa¬ 
tion of the generalized NFW profile has three solutions if 
1^1 < rji, where rj is the position of a source in the source 
plane and is a radius of the radial caustic in the source 
plane. The image separation is defined between the outer two 
solutions and is approximated as 9{M, zs, Zh) — 2^t/DoL, 
where is a radius of the tangential critical curve in the 
lens plane, Dql denotes the angular diameter distance from 
the observer to the lens plane, and zs and zl indicate the 
redshifts of the source and lens, respectively. 

The concentration parameter Cvir is one of the most 
important parameter in the generalized NFW density pro¬ 
file (eq. 0). In numerical simulations, it has been found 
that the concentration para meter depends on t he mass M 
and redshift z of halos fe.g.. lBullock et 11^ - Moreover, 
it shows considerable sc atter which reflects the difference 
in the formation epoch iWechsler et al l l2^ . The scatter 
in the concentration parameter is well described by a log¬ 
normal distribution. For the median of concentration pa¬ 
rameter Cvir .med, we adopt t he ma ss and redshift dependence 
reported bv iBullock et al.l (1200 Jl : 


Cvir.med (Tf, 2:) — (2 OL^- 


M 


1 -b 2 V 


( 2 ) 


Note that this fittin g form is slightly different from the one 
iBullock et al.l ^2001^ originally proposed, and is correct at 
M ~ Mq. a factor 2 — a in equation g ives a 

natural way to g eneralize a ^ 1 llKeeton fc Mad 3*111120011: 
l.ling fc Sutoll200al . The scatter of the concentration param¬ 
eter Oc is also important elem ent in gravitational lens statis¬ 
tics llKeeton fc Mada'nll2n0lh . We adopt the value ctc = 0.3 
has been obtained from N-body simu lations (Uinj 

Bullock et ani200ll: IWechsler et~al]|2002l: l.ling fc Sut3 

The probability that a source at redshift 2s and hav¬ 
ing the absolute luminosity L is observed as multiply lensed 
system with separation larger than 6 is given by 


which 

H ' 


p^i>e-zs,L) = 



dM 


^cdt ,, .3 dn 

^C'lensB-T (1 2 l) 

d2L 


dM’ 


where B is the magnification bias (E rne 


(3) 


y?^{zs,L) 



dyy^{zs,L/p{y)) 


1 


(4) 


with $(2s,I/) being the luminosity function of sources. The 
lensing cross section niens is simply given by the area encom¬ 
passed by the radial caustic, niens = '’^Vr Dqj^/Dq^, where 
Dos indicates the angular diameter distance from the ob¬ 
server to the source plane. The lower limit of mass integral 
Mmin is related to 0 as 0 zs, zl)- The magnifi¬ 

cation bias (eq. |3]) should be calculated for the faintest of 
the two images, because both lensed i mages must appea r 
above the flux limit of the 2QZ survey dMiller et alJl2f)oiill . 
We use an approxim ation of the magnif ication factor p{y) 
which was derived by lOguri et al.l (1200 2h . 

Since wide separation lenses with 0 > 30" are con¬ 
sidered to be generated by massive clusters, we should 
choose the mass function of dark halos ( dn/dM in eq. 
El) c arefully. We adopt equation (B3) of l.lenkins et all 
1I2OOIII whic h well agrees with simulated high-mas s halo 
abundan ce lEvrard et all 120021: iKomatsu fc S elia.kl I2OO2I : 
IHu fc KravtsovIPOO^^ier'Daoiret^Tl^Odll . Note that p = 
180fl(2)pcrit(z) should be used as the mean overdensity 
when one adopts this mass function. 


3 QUASAR CATALOG 


To make a precise prediction which can be directly com¬ 
pared with observed lensing rate, we must properly take 
account of the redshift and magnitude distributions. Since 
the whole sample used to search lensed quasars is not pub¬ 
licly available, we ins tead use the 2dF 10 k catalog compris¬ 
ing ~ 10000 quasars dCroom et alJl200ih which is a part of 
the whole sample. Then predicted numbers of wide separa¬ 
tion lensing are calculated as follows. First, from the cata¬ 
log we extract the numbers of quasars N{zi,bj) which are 
located Zi — Az/2 < Zi < Zi + Azl2 and have magnitude 
bj — Ab/2 < bj < bj + A&/2. We use A2 = 0.1 and Ab = 0.2. 
The average probability that quasars are lensed with sepa¬ 
rations larger than 0 is then given by 


Plens{>d) 




(5) 


where L{b) is the R-band absolute luminosity corresponds 
to b. When the quasar continuum spectrum is described by 
a power law, oc the K-correction can be approxi¬ 

mated as K{z) = —2.5(1—Qs) log(l-l-2). We assume Og = 0.5 
to calculate the K-correction. The total number of lensed 
quasars expected in the 2dF QSO survey is 


Nlens{>d) = NQSoPlena{>d), 


( 6 ) 


where Nqso = 22163 is the total number of quasars in the 
whole sample. 

The luminosity function of quasars is needed to com¬ 
pute magnification bias. We adopt the double power law 
luminosity function: 

( 2 , ^ [l/l4z)]9k + [l/l4z)]9i U{zy 
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Figure 1. The predicted number distribution of image separa¬ 
tions in the 2dF QSO survey (eq. |^). For the density profile of 
lens objects, both o = 1.5 and a = 1 are considered (see eq. 0). 
The numbers assuming Singular Isothermal Sphere (SIS) density 
profile, p oc are also plotted for reference. The cosmologi¬ 

cal models adopted in these plots are (OoiCrs) = (0.3, 0.7) and 
(Hq, (Tg) = (0.3,1.0) in the flat universe (Qq + -^o = !)• 


We assume pure luminosity evolution models with Lt,{z) oc 
-^Q kxz+k 2 z ^ ^ best-fitting model which was derived 

by iBovle et alJ J2nn(1i in the (no,Ao) = (0.3, 0.7) universe: 
Ph = 3.41, P, = 1.58, ki = 1.36, fcz = -0.27, and M. = 
-21.15 -b Slog/i. 

4 RESULTS 

4.1 Numbers of Lenses 

First we plot the predicted number distribution of image sep¬ 
aration in the 2dF QSO survey in Figure 0 The cases that 
the density profile of lens objects is described by the Singu¬ 
lar Isothermal Sphere (SIS) are also shown for reference. The 
value of as has been constrained from X-ray clusters or cos¬ 
mic shear, but resultant values show discrepancies among 
papers, ranging from ~ 0.7 to ~ 1.0 for Qq — 0.3 (e.g., 
IPiernaoli et alJl2003h . Therefore we plot both eg = 0.7 and 
(Tg = 1 models. Figure Q clearly indicates that the predicted 
numbers of lenses strongly depend on both density profile 
(a) and the abundance (erg) of dark halos. In particular, 
numbers of areminute-separation lensed quasars are highly 
sensitive to as because such very wide separation lenses are 
mainly produced by massive clusters. We find that expected 
numbers of lensed quasars with image separation larger than 
200” is quite small, A'’iens(> 200”) < 10“®, w hile 4 quasar 

S which are likely to be lensed were found iMiller et alJ 
. We further find that even an “extreme” model with 
a = 1.5 and as — 1.4 can produce only ~ 3 x 10“"^ lenses 
with 6 > 200” on average. On the other hand, less wide sep¬ 
aration lenses (9 ~ 30”) may be statistically possible if both 
a and erg are large. 

To see what constraint can be put from the detection 
of wide separation lenses, in Figure |21 we plot contours of 
iViens(> 30”) in the Qo-erg plane. We also plot constraints 
from the observed number of lenses assuming that some 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 

no=i-Ao 


Figure 2. Contours of N\ens{> 30”) in the flo-crs plane. Con¬ 
straints (95% confidence level) from the detection of wide sepa¬ 
ration lenses are shown by dashed lines. From lower to upper of 
dashed lines, the number of genuine lens systems is assumed to 
be 1, 2, 3, 4, 5, 6. Dotted lines indicate the recent constraints 
on Qq and eg from X-ray clusters or cosmic shear which can be 
approximately expressed as 0.4 < o-gfl^ ® < 0.55. 

of six pairs reported by iMiller et alJ ^2003^ are true lens 
systems. These constraints are calculated from the Poisson 
distribution P{k\N) = N^e~^/k\, where k is the observed 
number and N is the expectation. Given the numbers of 
true lenses range from fc = 1 to 6, the lower limits of N 
(95% confidence level) become 0.35, 0.82, 1.36, 1.97, 2.61, 
and 3.29, respectively. This figure indicates that the detec¬ 
tion of even one lens with 9 > 30” in the 2dF QSO catalog 
needs very high as, as > 1.2 for flo = 0.3 and a = 1.5. 
Our constraint and constraints from X-ray clusters or cosmic 
she ar which are approxim ated as 0.4 < crgQo ® < 0.55 (see 
e.g.. lPierDaoli et al.l2003|l shows marginal agreement if Qo ^ 
0.15, as > 1.5, and a = 1.5, although this solution seems 
too extreme. In particular, models with a = 1.0 are hardly 
acceptable because they need unusually high erg or Qq in or¬ 
der to produce lenses with 9 > 30”. However, this result is 
somewhat embarrassing bec ause the lack of wide separation 
lensing in other surveys (e.g.. lPhilliDS. Browne, fc WilkinsonI 
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Table 1. Theoretical uncertainties in the expected numbers of 
lensed quasars. The fiducial model has a. = 1.5 and (Qo 5 <^ 8 ) = 
(0.3,1.0). “Press Schechter MF” and “Sheth Tormen 
MF” models are same as the fiducial model exc e pt fo r us- 
inff the mass fun c tion d erived by IPress fc Schecht'e^ _l|lj}7jj_and 
ISh gUi &: TormeJ instead of that fitted hT TjeiildTisei~an 

J200l|), respectively. In “Groom et al. LF” modelj_we_use_tte 

S ir luminosity function from the 10 k catalog |^room_et_aJj 
. Note that this luminosity function has fairly shallower 
bright- and faint-end slopes, = 3.28 and 0i = 1.08. 


Model 

Nlen.(> 30") 

200 ") 

Fiducial Model 

7.9 X 10-2 

1.8 X lO-l" 

Press & Schechter MF 

1.2 X 10-1 

1.6 X lO-l" 

Sheth & Tormen MF 

2.8 X 10-1 

8.3 X 10-® 

Croom et al. LF 

4.4 X 10-2 

1.1 X lO-l" 


l2001al : IPhillips et alJl2001bl : lOfek et al.lf200li . l2003l already 
puts the upper limit of the lensing rate. For instance, our 
model with a — 1.5, flo = 0.3, and as = 1.2 predicts 
a lensing rate P{> 6”) ~ 3 x 10~^ at 2 ~ 1.3 which is 
margin ally consistent with t he upper limit of the lensing 
rate in IPhillips et a 0 (I2OOIIJI . P{> 6") < 3 X 10"'^ at 95% 
confidence limit. 

Note that these constraints are for the case that only 
one of six quasar pairs is a genuine lens system; if the num¬ 
ber of lenses is more than one, the situation becomes worse. 
In this case, the discrepancy between strong lensing con¬ 
straints and X-ray/shear constraints becomes more serious. 
Moreover, this result may conflict with other surveys which 
could not detect wide separation lenses. 

4.2 Theoretical Uncertainties 

We also examine possible theoretical uncertainties except for 
the uncertainties of the density proHle and cosmological pa¬ 
rameters. More specifically, we examine the uncertainties of 
the halo mass function and the quasar luminosity function. 
The result i s summarized in TableU F irst we adopt the mass 
functi ons of IPress fc Schechteil l^l^4^ an d ISheth fc Torm*^ 
^1999^ and see how the number of lensed quasars changes. 
We find that the uncertainty of the mass function is fairly 
large. In particular, the mass function of ISheth fc TormerJ 
^1999^ predicts more than three times as large number of 
lenses with 6 > 30" as our fiducial model. This is because 
the mass function of ISheth fc TormerJ ^1999^ seems to over¬ 
estimate the number densit y of massive halos il.Tenkins et al.l 
|2001I : IHu fc Kravtsovl20‘n^ . Moreover, we point out that the 
number of anomalously wide separation lenses {6 > 200") 
is much more sensitive to the choice of the mass function. 
Therefore, in the statistics of such anomalously wide sepa¬ 
ration lenses the uncertainty of the mass function should be 
carefully examined. 

Since the magnification bias is sensitive to the slope of 
adopted quasar luminosity function, next we examine the 
uncertainty of the quasar luminosity fu nction using the lu - 
minosity function from the 10k catalog dCroorn et alJlioOlll . 
This luminosity function has somewhat shallower slopes, 
Ph = 3.28 and j3i = 1.08, compared with our fiducial model. 
We find that this uncertainty is less than factor 2 and is not 
so large as to change our main results, because a, flo, and as 



1 10 ' 10 ^ 10 " 
At[yr] 


Figure 3. The conditional probability distributions of possible 
time delays are shown for each quasar pain_The_^finiton and 
calculation of P{> At\9) was shown by |Oguri_et_^ J200^. The 
redshift and magnitude for each image are taken into account. 
For the cosmological model, we adopt {Qf),crs) = (0.3,1.0). 

can change the number of wide separation lenses by orders 
of magnitude. 

4.3 Expected Time Delays between Images 

The main drawback of wide separation lensing statistics is 
that it is hard to recognize quasar pairs as gravitational lens 
systems due to large time delays between images; any spec¬ 
trum change within the timescale of differential time delays 
may prevent one from selecting such systems as lens candi¬ 
dates. We show con ditional probability distributions of time 
delays proposed bv lOguri et alJ (120021) in Figure |3| We find 
that the dependence of time delays on the density profile is 
weak when the separat ion is large, althou gh this tendency 
was already shown by lOguri et alJ l^2002^ . They also con¬ 
cluded that time delay probability distribution is insensitive 
to cosmological parameters. Figure suggests that lenses 
with even smaller separations, 6 ~ 30", are likely to have 
time delays larger than ten years which may be typ ical time 
scale forming the broad absorption line (lMall200^ . There¬ 
fore, to assert quasar pairs as lensed images, one must com¬ 
pare spectral signature which would be unchanged within 
possible time delays. Needed time scales are At > lOOyr for 
9 ~ 30" and At > lOOOyr for 9 ^ 200". We note that the 
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information of flux ratio and the central core imag e may 
becom e one of the evidence of gravitational leasing jRusinl 

li^. 


5 SUMMARY 

In this paper, we have calculated the predicted numbers 
of arcminute-separation {9 > 30”) lensed quasars in the 
2dF QSO survey. We have presented realistic predictions 
based on the CDM scenario taking account of the redshift 
and magnitude distributions of the quasar catalog. Detailed 
comparison between theoretical and observed numbers of 
lensed quasars indicates that the detection of wide separa¬ 
tion lenses puts interesting constraints on the density prohle 
and abundance of dark halos. The case that only one of six 
pairs is genuine lens system is marginally consistent with 
the model that has cuspy inner density profile p oc r~^'^ 
and the large value of as, ag > 1.2 for flo = 0.3. To rec¬ 
oncile this result with X-ray or shear measurement, much 
smaller flo (Ho ^ 0.15) and much larger ag (ag > 1.5) 
are needed. Our result of this large ag is similar to that 
of the Sunyaev-Zekdovich angular power spectrum, ag ~ 1.1 
iKomatsu fc Selia 31^. We have found also that it is quite 
hard to produce lenses with separation > 200”. Thus a con¬ 
servative interpretation of this observation is that none of 
these quasar pairs is lensed. But if it turns out that some of 
these quasar pairs are genuine lens systems, we can put in¬ 
teresting constraints on not only the density profile of dark 
halos but also the Ho and ag that are somewhat different 
from X-ray/shear constraints. We note that the number of 
genui ne lens systems ma y be significantly larger than six 
which iMiller et al.l (l2003ll reported because only 11 of 38 
candidates has been observed spectroscopically. However, 
such anomalously high leasing rate cannot be reproduced 
by even the most optimistic models. In this case, we have to 
examine whether we miss some important systematic effects 
which increase lensing rates. One possible systematic effect 
is the asymmetry of lensing halos, although this effect has 
been considered to be small so far. In any case, statistics of 
wide separation lensing offer a promising way to probe the 
abundance and density profile of dark halos. 
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